Unified models of the cosmological dark sector 
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Abstract. We model the cosmological substratum by a viscous fluid that is supposed to 
provide a unified description of the dark sector and pressureless baryonic matter. In the 
homogeneous and isotropic background the total energy density of this mixture behaves as 
a generalized Chaplygin gas. The perturbations of this energy density are intrinsically non- 
adiabatic and source relative entropy perturbations. The resulting baryonic matter power 
spectrum is shown to be compatible with the 2dFGRS and SDSS (DR7) data. A joint statistical 
analysis, using also Hubble-function and supernovae la data, shows that, different from other 
studies, there exists a maximum in the probability distribution for a negative present value 
go ~ —0.53 of the deceleration parameter. Moreover, different from other approaches, the 
unified model presented here favors a matter content that is of the order of the baryonic matter 
abundance suggested by big-bang nucleosynthesis. 



1. Introduction 

According to the prevailing interpretation, our Universe is dynamically dominated by a 
cosmological constant A (or a dynamical equivalent, called dark energy (DE)) which contributes 
more than 70% to the total cosmic energy budget. More than 20% are contributed by cold dark 
matter (CDM) and only about 5% are in the form of conventional, baryonic matter. Because 
of the cosmological constant problem in its different facets, including the coincidence problem, 
a great deal of work was devoted to alternative approaches in which a similar dynamics as that 
of the ACDM model is reproduced with a time varying cosmological term, i.e., the cosmological 
constant is replaced by a dynamical quantity. Both dark matter (DM) and DE manifest 
themselves so far only through their gravitational interaction. This provides a motivation for 
approaches in which DM and DE appear as different manifestation of one single dark-sector 
component. The Chaplygin-gas model and its different generalizations realize this idea. Unified 
models of the dark sector of this type are attractive since one and the same component behaves 
as pressureless matter at high redshifts and as a cosmological constant in the long time limit. 
While the homogeneous and isotropic background dynamics for the (generalized) Chaplygin 
gas (GCG) is well compatible with the data, the study of the perturbation dynamics resulted 
in problems which apparently ruled out all Chaplygin-gas type models except those that are 
observationally almost indistinguishable from the ACDM model. To circumvent this problem, 
nonadiabatic perturbations were postulated and designed in a way to make the effective sound 
speed vanish. But this amounts to an ad hoc procedure which leaves open the physical origin 
of nonadiabatic perturbations. There exists, however, a different type of unified models of 



the dark sector, namely viscous cosmological models, which are intrinsically nonadiabatic [T]. 
In the homogeneous and isotropic background a one-component viscous fluid shares the same 
dynamics as a GCG. Now, what is observed in the redshift surveys is not the spectrum of the 
dark-matter distribution but the baryonic matter spectrum. Including a baryon component into 
the perturbation dynamics for a universe with a Chaplygin-gas dark sector, there appears the 
new problem that the unified Chaplygin-gas scenario itself is disfavored by the data. It is only 
if the unified scenario with a fixed pressureless (supposedly) baryonic matter fraction of about 
0.043 (according to the results from WMAP and primordial nucleosynthesis) is imposed on the 
dynamics, that consistency with the data is obtained [2]. If the pressureless matter fraction is 
left free, its best-fit value is much larger than the baryonic fraction. In fact it becomes even 
close to unity, leaving only a small percentage for the Chaplygin gas, thus invalidating the entire 
scenario. In other words, a Chaplygin-gas-based unified model of the dark sector is difficult to 
reconcile with observations. One may ask now, whether the status of unified models can again be 
remedied by replacing the Chaplygin gas by a viscous fiuid. It is exactly this question that we are 
going to investigate in the present paper. It is our purpose to study cosmological perturbations 
for a two-component model of baryons and a viscous fiuid, where the latter represents a one- 
component description of the dark sector [3]. We shall show that such type of unified model 
is not only consistent for a fixed fraction of the baryons but also for the case that the matter 
fraction is left free. Our analysis demonstrates that the statistically preferred value for the 
abundance of pressureless matter is compatible with the mentioned baryon fraction 0.043 that 
follows from the synthesis of light elements. 



2. The two-component model 

The cosmic medium is assumed to be describable by an energy-momentum tensor T*'^ which 
splits into a matter part T}^ and viscous fiuid part Ty, 

= pu'v!" +p(g''' + u'u^) , T^ = T^ + T^, (1) 

with 

T^i = Pmv^mUm + PM (g'^ + uiiult) , = pvv^.vl'y + pv (g'^ + u\rv!f) , (2) 

where the subscript "M" stands for matter and the subscript "V" stands for viscous. The total 
cosmic fluid is characterized by a four velocity while Uj^^ represents the four velocity of 
the matter part and Uy represents the four velocity of the viscous fluid. Energy- momentum 
conservation is supposed to hold separately for each of the components, 

t!}., = T{,'., = ^ T% = 0. (3) 

Let the matter be pressureless, i.e. pM = and the viscous fluid be characterized by a bulk 
viscous pressure pv = p = — C©) where C = const and = n*^ is the fluid expansion. Under this 
condition the total pressure coincides with the pressure of the viscous component. In terms of 
the present value qo of the deceleration parameter the Hubble rate can be written as [T] 

^ ^ 1 - 2(70 + 2 (1 + go) , (4) 
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which coincides with the Hubble rate of a specific (a = —1/2) GCG. Since pM = PMoa~'^, we have 
Pv = p — pmoo,~^ . It is the total energy density that behaves as a GCG, not the component V . 
This type of unified model differs from unified models in which the total energy density is the sum 
of a GCG and a baryon component. Only if the baryon component is ignored, both descriptions 



coincide. Consequently, in the homogeneous and isotropic background, a generahzed Chaplygin 
gas with a = —1/2 can be seen as a unified description of the cosmic medium, consisting of 
a separately conserved matter component and a bulk viscous fluid with = const, where the 
latter itself represents a unified model of the dark sector. 

3. Perturbations 

The nonadiabaticity of the system as a whole is characterized by 

^ ^ p m^i^-^ . (5) 



P + P PP^V p \p 

The quantity (0) is governed by the dynamics of the total energy-density perturbation p and 
by the perturbations of the expansion scalar, which is also a quantity that characterizes the 
system as a whole. The behavior of these quantities is described by the energy-momentum 
conservation for the entire system and by the Raychaudhuri equation, respectively. Both of 
these equations are coupled to each other. The remarkable point is that these quantities and, 
consequently, the total energy density perturbation, are independent of the two-component 
structure of the medium. The reason is the direct relation p = —(^0 between the pressure 
perturbations and the perturbations of the expansion scalar. It is convenient to describe the 
perturbation dynamics in terms of gauge invariant quantities which represent perturbations on 
comoving hyper surf aces, indicated by a superscript c. These are defined as (v is the velocity 
potential, defined by u^ = 

P^ _ P s: _ /''^ 



Then, in linear order, 



+ v, 6 = ^. (6) 
P P P 

5" + f{a)6' + g{a)5 = Q, (7) 

where 5' = ^ and the coefficients / (a) and g (a) are given by explicitly known background 
quantities [UE]. The relative entropy perturbations are defined by 

Smv^'-^-^^ (8) 

PM PV + PV 



and obey the equation 



S'vM + r{a)SvM + s{a)SvM = c{a)6' + d{a)5 . (9) 



Again, the coefficients r{a), s{a), c(a) and d{a) are combinations of analytically given 
background variables [3j. The quantity relevant for the observations is the fractional 



perturbation 5m = ^ of the energy density of the baryons, given by 



Sm = - 

7 



jr PV+Pq 

bvM 

P 



(10) 



At early times, i.e. for small scale factors a <C 1, the equation ([7]) has the asymptotic form 
5" + ^ (5' — 2^2- (5 = independent of go and for all scales. Eq. ([9]) then reduces to 'S'v'm + ^'^V'M ~ 
with the solution SyM = const = 0. Consequently, there are neither nonadiabatic contributions 
to the total energy-density fiuctuations nor relative entropy perturbations and we have purely 
adiabatic perturbations 5m = (5 at a <C 1. This allows us to match the initial conditions for 
equations d?]) and dH) with those of the ACDM model at early times. 
The results of our statistical analysis are shown in Fig. [TJ 




Figure 1. Left panel: Hubble diagram, center panel: probability distribution function (PDF) for qo, based on 
a joint analysis of matter power spectrum and SN la data. The right panel shows the PDF for the pressureless 
component. For more details see [3]. 

4. Discussion and Conclusions 

The approach presented here is based on the fact that the two-component system of a bulk 
viscous fluid and a separately conserved baryon component behaves in the background as a 
generalized Chaplygin gas with a = While the baryon component may be considered 

dynamically negligible in the background, the situation is different on the perturbative 
level, since the observed matter agglomerations are related to baryonic density fluctuations. 
These fluctuations are obtained from a combination of the nonadiabatic total energy density 
perturbations and relative entropy perturbations in the two-component system where the former 
source the latter. The probability distribution for the deceleration parameter has a maximum 
at qo ~ —0.53 which partially removes the degeneracy of previous studies which, taken at face 
value, were incompatible with an accelerated expansion and thus in obvious tension with results 
for the background. Perhaps still more important is the test of the unified model itself. Many 
investigations on approaches with a unified dark sector fix the pressureless matter component 
to be that of the favored (by the WMAP data) baryon fraction and then check whether or not 
the resulting dynamics can reproduce the observations. But this does not say anything on how 
probable the division of the total cosmic substratum into roughly 96% of a dark substance and 
roughly 4% of pressureless matter is. To decide this question, one has to consider the pressureless 
matter fraction as a free parameter and to find out which abundance is actually favored by the 
data. Our analysis revealed that the matter fraction probability is indeed highest for values 
smaller than roughly 8%. This is a result in favor of the unified viscous model. We recall that 
a corresponding analysis for a Chaplygin gas results in values close to unity [2J which seems to 
rule out such type of approaches. The present viscous model, on the other hand, remains an 
option for a unified description of the dark sector, at least as far as the matter power spectrum 
is concerned. 
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